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Abstract 

i <~| A remarkable difference between quantum and classical programs is that the con- 

Q-i trol flow of the former can be either classical or quantum. One of the key issues in 

the theory of quantum programming languages is defining and understanding quantum 
control flow. A functional language with quantum control flow was defined by Al- 
tenkirch and Grattage [Proc. LICS'05, pp. 249-258]. This paper extends their work, 
and we introduce a general quantum control structure by defining three new quantum 
program constructs, namely quantum guarded command, quantum choice and quantum 
recursion. We clarify the relation between quantum choices and probabilistic choices. 
An interesting difference between quantum recursions with classical control flows and 
with quantum control flows is revealed. 

1 Introduction 

Since Knill [ 8 ] introduced the Quantum Random Access Machine (QRAM) model for quan- 
tum computing and proposed a set of conventions for writing quantum pseud-ocodes in 
1996, several quantum programming languages have been defined in the last 16 years; for 
example QCL by Omer [12J, qGCL by Sanders and Zuliani [ 13 ], QPL by Selinger lfl4l . and 
see Q for an excellent survey. One of the key design ideas of almost all existing quantum 
languages can be summarised by the influential slogan "quantum data, classical control" 
proposed by Selinger lfl4ll . meaning that the control flow of a quantum program is still clas- 
sical, but the program operates on quantum data. An exception is Altenkirch and Grattage's 
functional language QML [2], where "quantum control" flow was introduced; more pre- 
cisely, they observed that in the quantum setting the case construct naturally splits into two 
variants: 

• case, which measures a qubit in the data it analyses; 

• case , which analyses quantum data without measuring. 

The control flow in the case construct is determined by the outcome of a measurement and 
thus is classical. However, a quantum control flow appears in the case construct, as shown 
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in the following example where a special form of the case , namely the if° — then — else 
statement, is used. 

Example 1.1 Basic quantum gates implemented in QML [2]: The Hadamard gate is writ- 
ten as: 

had : Q2 — ° Q2 
had x = if x 

then {— ^=(qf alse — qtrue)} 

else {— p(qf alse + qtrue)} 

V2 

and the CNOT gate is as follows: 

cnot : Q 2 — ° Q2 — Q2 <8> Q2 
cnot c x = if c 

then (qtrue, not x) 

else (qfalse,x) 

where Q2 is the type ofqubits, and not is the NOT gate: 

not : Q 2 — o Q 2 
not x = if x 

then qf alse 

else qtrue 

A new research line of quantum programming with quantum control flow was then 
initiated by Altenkirch and Grattage in and further pursued by themselves and others in 
a series of papers J3][9]|. 

The present paper continues this line of research, and we extend the idea of "quantum 
control" by introducing three new quantum program constructs: 

(1) Quantum Guarded Command: Our first step toward a general quantum control 
structure is to introduce a quantum generalisation of Dijkstra's guarded command [6]. Re- 
call that a guarded command can be written as follows: 

□?=i h -> Ct (i) 

where for each 1 < i < n, the command C{ is guarded by the Boolean expression b{, and Cj 
will be executed only when b{ is true. Obviously, the case operator in QML is a quantum 
generalisation of guarded command with classical control. On the other hand, as shown 
in the above example, the case operator in QML implements a unitary transformation by 
decomposing it into two orthogonal branches along the quantum control flow determined 
by a chosen qubit. So, it is already a kind of guarded command with quantum control flow. 

An even clearer idea for defining quantum guarded command stems from a quite differ- 
ent area, namely quantum walks |4ll, ifTI: 
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Example 1.2 Quantum walks on graphs /[7]/. - Let (V, E) be an n— regular directed graph. 
Then we can label each edge with a number between 1 and n such that for each 1 <i <n, 
the directed edges labeled iform a permutation. Let T~Ly be the Hilbert space spanned by 
states {\v)} V £.y. Then for each 1 < i < n, we can define a shift operator Si on T~Ly: 

Si\v) = | the ith neighbour of v) 

for any v G V. We introduce an auxiliary quantum variable q with the state Hilbert space 
% q spanned by {|i)}" =1 . Now we are able to combine the operators Si (1 < i < n) along 
q to form a whole shift operator: 

S = U r U q, \i) -> Si (2) 

on Ti q <g> Tiy: 

S\v,i) = (Si\v))\i) (3) 

for any 1 < i < n and v £ V. If we further choose a unitary operator U on T-L q then a 
coined quantum walk on graph (V, E) is defined by modelling its single step by the unitary 
operator: 

W = S{I U y ® U) 

on Tiy (8) T~L q , where I% v is the identity operator in Tiy. Usually, H q is called the "coin 
space", and U the "coin-tossing operator". 

The guarded command notation is adopted in Eq. (O to indicate that the shift opera- 
tor S is indeed a guarded command with quantum control. It is interesting to note that 
both Examples 11.11 and 11.21 defined a guarded command with quantum control, but their 
defining strategies are quite different: in Example 11.11 a quantum control flow is detected 
by decomposing a unitary operator along an existing qubit; in contrast, a quantum control 
flow is created in Example 1 1.2l by introducing a new quantum variable so that we can com- 
bining a family of unitary operators along the created flow. The defining strategy used in 
Example 1 1 . 21 naturally leads us to a general form of quantum guarded command: 

□tl q, \i) -> Pi (4) 

where Pi, P n are a family of quantum programs, and a new family of quantum variables 
q that do not appear in Pi , . . . , P n is introduced so that we can form a quantum guarded 
command by combining P\, ...,P n along an orthonormal basis {\i)} of the state space of q. 
For each 1 < i < n, Pi is guarded by the basis state and a superposition of these basis 
states yields a quantum control flow. 

(2) Quantum Choice: Guarded commands are the most widely accepted mechanism 
for nondeterministic programming. Nondeterminism in guarded command £T|) is a conse- 
quence of the "overlapping" of the guards b\,...,b n . In particular, if b\ = ■ ■ ■ = b n = true, 
then guarded command (fl]) becomes a demonic choice: 

□? =1 C h (5) 
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where the alternatives are chosen unpredictably. Usually, the demonic choice is sepa- 
rately defined as an explicit operator rather than a special case of guarded command due to 
its importance as a means of abstraction in programming. 

To formalise randomised algorithms, research on probabilistic programming [10] started 
in 1980's with the introduction of probabilistic choice: 

□r=i Ci@pi, (6) 

where {p^ is a probability distribution; that is, Pi > for all i, and Y^l=i Pi = !• The 
probabilistic choice © randomly chooses the command C, with probability pi for every i, 
and thus it can be seen as a refinement of the demonic choice ©. A probabilistic choice is 
often used to represent a decision in forks according to a certain probability distribution in 
a randomised algorithm. 

A natural question then arises in the realm of quantum programming: is it possible to 
define a quantum choice of programs? Indeed, an idea is already there in the construction 
of quantum walks, although not explicitly stated. In Example 11.21 each shift operator Si 
can be considered as an independent program, the "coin-tossing operator" U is employed to 
create a superposition of Si (1 < i < n), and thus the single step operator W can be seen as 
a quantum choice among Si (1 < i < n). Extending the idea used in Example 1 1.21 we can 
define a general quantum choice as a sequential composition of a "coin-tossing" program 
and a quantum guarded command: 

□? =1 P;q, \i) -> P = P; D" = i q, \i) -> Pi, (7) 

where Pi, P n are a family of quantum programs, q is a new family of quantum variables 
with {\i}} as an orthonormal basis of its state space, and P is a quantum program acting 
on q. Intuitively, quantum choice ((TJ) first runs program P to produce a superposition of the 
execution paths of programs Pj (1 < i < n), and then the guarded command D™ =1 q,\i) —> 
Pi follows. During the execution of the guarded command, each Pj is running along its 
own path within the whole superposition of execution paths of Pj (1 < i < n). It is 
widely accepted that quantum superposition is responsible for the advantage of quantum 
computers over classical computers. The power of superposition of quantum states has been 
successfully exploited in quantum computing. Quantum choices may provide a platform for 
explore a higher level of quantum superposition in computing, namely the superposition of 
quantum programs. 

(3) Quantum Recursion: Most classical programming languages allow direct specifi- 
cation of recursive procedures. Quantum loops and more general quantum recursive proce- 
dures were already defined in Selinger's language QPL lfl4ll . and termination of quantum 
loops were analysed by the authors in lfT8l . But quantum recursions considered in lfl4l[T8Ti 
contain no quantum control flows because there branchings in quantum programs are all de- 
termined by the outcomes of quantum measurements. After introducing quantum guarded 
commands and quantum choices, loops and recursive procedures with quantum control 
flows can be defined. As will be seen later, a major difference between quantum recursions 
with and without quantum controls is: auxiliary quantum variables must be introduced in or- 
der to define quantum controls. Thus, localisation mechanism is needed in defining quantum 
recursions with quantum control so that consistency of quantum variables is guaranteed. 
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1.1 Technical Contributions of the Paper 

As shown above, a general notion of quantum control flow comes naturally out from gen- 
eralising the case construct in Altenkirch and Grattage's language QML and the shift 
operators in quantum walks. However, a major difficulty arises in defining the semantics of 
quantum guarded commands. For the case where no quantum measurement occur in any Pj 
(1 < i < n), the semantics of each P{ is simply a sequence of unitary operators, and the se- 
mantics of guarded command © can be defined in exactly the same way as Eq. (0. When- 
ever some Pi contains quantum measurements, however, its semantic structure becomes a 
tree of linear operators with branching happening at the points where measurements are 
performed. Then defining the semantics of guarded command ((U) requires to properly com- 
bine a collection of trees such that certain quantum mechanical principles are obeyed. This 
problem will be circumvented in Sec. [3j 

1.2 Organisation of the Paper 

A new quantum programming language QGCL with quantum guarded commands is defined 
in Sec. [2 Sec. [Jjprepares some key ingredients needed in defining the semantics of QGCL. 
The denotational semantics and weakest precondition semantics of QGCL are presented 
in Sec. [4] In Sec. [5j quantum choice is defined in terms of quantum guarded command, 
and probabilistic choice is implemented by quantum choice by introducing local variables. 
Because of the limited space, quantum recursion is only briefly touched in Sec. [6l For 
readability, all proofs are deferred to the Appendix. 

2 QGCL: A Language with Quantum Guarded Commands 

We now define a quantum programming language QGCL with quantum guarded commands. 
QGCL is essentially an extension of Sanders and Zuliani's qGCL obtained by adding quan- 
tum control flow. But the presentation of QCGL is quite different from qGCL due to the 
complications in the semantics of quantum guarded commands. We assume a countable set 
qVar of quantum variables ranged over by q, qi, q<z, .... For simplicity of the presentation, 
we only consider a purely quantum programming language, but we include a countably in- 
finite set Var of classical variables ranged over by x, y, ... so that we can use them to record 
outcomes of quantum measurements. However, classical computation described by, for ex- 
ample, the assignment statement x := e in a classical programming language is excluded. 
It is required that the sets of classical and quantum variables are disjoint. For each classical 
variable x £ Var, its type is assumed to be a non-empty set D x ; that is, x takes values from 
D x . For each quantum variable q G qVar, its type is a Hilbert space type(q) = 7-L q , which 
is the state space of the quantum system denoted by q. For a sequence q = q\ , q2, ■ ■ ■ of 
quantum variables, we write: 

ty P e(q)=H4 = (g>H gi . 

i>l 
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Similarly, for any set V C qVar, we write: 

type{V)=H v = 0W g . 

In particular, we write H a « for type(qVar). To simplify the notation, we often identify a 
sequence of variables with the set of these variables provided they are distinct. 

Definition 2.1 For each QGCL program P, we write var(P) for the set of its classical 
variables and qvar(P) for its quantum variables. QGCL programs are inductively defined 
as follows: 

1. abort and skip are programs, and 

var (abort) = var(skip) = 0, 
gvar(abort) = qvar(skip) = 0. 

2. Ifq = qi,...,qkisa sequence of quantum variables, and U is a unitary operator on 
type(q), then U[q] is a program, and 

var(U\q\) = 0, qvar(U[q}) = q. 

3. Ifq = q\,...,qkisa sequence of quantum variables, x is a classical variable, M = 
{M m } is a quantum measurement in type(q), and {P m } is a family of programs 
indexed by the outcomes m of measurement M such that spec(M) C D x , where 
spec(M) = {m} is the spectrum of M; that is, the set of all possible outcomes of M, 
and x ^ |J m var(P m ), then 

P = M[x <- q] : {P m } (8) 

is a program, and 

var(P) = {x} U\J var (P m ), 

m 

qvar(P) = q U qvar{P m ). 

m 

4. Ifq = qi, (j% is a sequence of quantum variables, {|i)}™ =1 is an orthonormal basis 
oftype(q), and {Pi}f =1 is a family of programs such that 

n 

qn[j qVar{Pi) = 0, 
i=i 

then 

P = \i)^Pi 
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is a program, and 

n 

var(P) = var(Pi), 
i=i 

n 

qvar(P) = qU qvar(Pi). 

i=l 

5. If P\ and P2 are programs such that var{P\) n variP-z) = 0, then P\; P2 is a pro- 
gram, and 

var(P\\ P2) = var(Pi) U var(P2), 
qvar(Pi; P2) = qvar{P\) U qvar(P2). 

The meanings of abort and skip are the same as in a classical programming language. 
Two kinds of statements are introduced in the above definition to describe basic quantum 
operations, namely unitary transformation and measurement. In the unitary transformation 
U[q], only quantum variables q but no classical variables appear, and the transformation is 
applied to q. In statement ([8]), a measurement M is first performed on quantum variables q 
with the outcome stored in classical variable x, and then whenever outcome m is reported, 
the corresponding subprogram P m is executed. The intuitive meaning of quantum guarded 
command was already carefully explained in Sec. [T] Whenever the sequence q of quantum 
variables can be recognised from the context, Df =1 q, \i) — > Pi can be abbreviated to 
□f =1 \i) —7- P{. The sequential composition Pi; P2 is similar to that in a classical language, 
and the requirement var(Pi) n var(P2) = means that the outcomes of measurements 
performed at different points are stored in different classical variables. Such a requirement 
is mainly for technical convenience, and it will considerably simplify the presentation. The 
syntax of QGCL can be summarised as follows: 

P := abort | skip | P\] P2 

\U\q] (unitary transformation) 

I measure M[g] : {Pm} (quantum measurement (9) 

= classical guarded command) 
I □"=! q,\i) — > Pi (quantum guarded command) 

3 Guarded Compositions of Quantum Operations 

3.1 Guarded composition of unitary operators 

A major difficulty in defining the semantics of QGCL comes from the treatment of guarded 
commands where a guarded composition of semantic functions is vital. To ease the under- 
standing of a general definition of such a guarded composition, we start with the guarded 
composition of unitary operators, which is a straightforward generalisation of the quantum 
walk shift operator S in Example 1 1.21 
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Definition 3.1 For each 1 < i < n, let Ui be an unitary operator in Hilbert space %. Let 
% s be a Hilbert space with {\i)} as an orthonormal basis. Then we define a linear operator: 



i=l 



in %<®1-L S by 



u{\m)) = (m))\i) 

for any G H and for any 1 < i < n. Then by linearity we have: 

/ n \ n 



(10) 



vi=l 



i=l 



for any \ipi), \ip n ) G %■ The operator U is called the guarded composition ofUi (1 < 
i <n) along {\i)}. 

Example 3.1 Quantum multiplexor: As a quantum generalisation of multiplexor, a well- 
known notion in digit logic, quantum multiplexor ( QMUXfor short) was introduced in U5\l 
as a useful tool in synthesis of quantum circuits. A QMUX U with k select qubits and 
d—qubit-wide data bus can be represented by a block-diagonal matrix: 



( U 



\ 



Ui 



U = diag(U , U\, U 2 k_ 1 ) 
Multiplexing Uo, U\, U 2 k_ 1 with k select quits is exactly the guarded composition 

along the computational basis of k qubits. 



Lemma 3.1 The guarded composition Df =1 \i) — > Ui is an unitary operator in 1-1®% 



3.2 Operator-valued functions 

For any Hilbert space H, we write for the space of (linear) operators on H. 

Definition 3.2 Let A be a nonempty set. Then a function F : A — )■ C{%) is called an 
operator-valued function in % over T, if 

^2F{5^-F(S)ni u , (11) 

<5eA 

where 1% is the identity operator in %, and C stands for the Lowner order; that is, A C B 
if and only if B — A is a positive operator. In particular, F is said to be full when Eq. f li 71 ) 
becomes equality. 
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The simplest examples of operator-valued function are unitary operators and measure- 
ments. 

Example 3.2 1. A unitary operator on Hilbert space % can be seen as a full operator- 
valued function over a singleton A = {e}. 

2. A measurement M on Hilbert space % can be seen as a full operator-valued function 
over its spectrum Spec(M). 

More generally, a super-operator defines a family of operator-valued functions. Let £ 
be a super-operator on Hilbert space %. Then £ has the Kraus operator-sum representation: 
£ = Ei o e\, meaning: £(p) = ^ EipEj for all density operators p 'mJi. For such 
a representation, we set A = {%] for the set of indexes, and define an operator-valued 
function over A by F(i) = Ei for every i. Since operator-sum representation of £ is 
not unique, £ defines not only a single operator-valued function. We write ¥{£ ) for the 
family of operator-valued functions defined by all Kraus operator-sum representations of £. 
Conversely, an operator-valued function determines uniquely a super-operator. 

Definition 3.3 Let F be an operator-valued function in Hilbert space % over set A. Then 
F defines a super-operator £{F) in % as follows: 

= J>(*)oF(*)t. 

For a family F of operator- valued functions, we write £{¥) = {£(F) : F G F}. It is 
obvious that £(¥(£)) = {£}. On the other hand, for any operator-valued function F over 
A = {<5i, ...,5^}, Theorem 8.2 in [ 11 ] indicates that ¥(£(F)) consists of all operator- valued 
functions G over T = {71, 7;} such that 

n 

G{ ll ) = Y,u l3 -F{5 3 ) 
i=i 

for each 1 < i < n, where n = max(fc, I), U = (uij) is an n x n unitary matrix, F(5i) = 
G(jj) = Oy, for all k + 1 < i < n and I + 1 < j < n. 

3.3 Guarded composition of operator- valued functions 

We first introduce a notation. Let Aj be a nonempty set for every 1 < i < n. Then the 
superposition of Aj (1 < i < n) is defined as follows: 

n 

Aj = {0™ = i<5j : Si e Aj for every 1 < i < n}. 

i=l 

Definition 3.4 For each 1 < i < n, let Fi be an operator-valued function in Hilbert space 
T~L over set Aj. Let iH s be a Hilbert space with {\i)} as an orthonormal basis. Then the 
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guarded composition of Fi (1 < i < n) along {\i)} is defined to be the operator-valued 
function in % (g) % s over ®" =1 Aj/ 

f =\ nu \i) F it 

F(®tA) (ei^>k>) = E (n A ^ ) (mmw) on 

\i=l / i=l \k^i J 

for any \ip n ) S % and for any 5i G Aj (1 < i < n), where 



/« particular, if F k is full and d = dim% < oo, 



, _ trF k {5 k ^F k (5 k ) 

" \H Tk ^ k trF k (r k fF k (r k y ^ 



_ t trF k {5 k fF k {5 k 



d 

for any 5 k £ A k (1 < k < n). 

It is easy to see that whenever Aj is a singleton for all 1 < i < n, then Eq. (fl2l 
degenerates to Eq. (fTOl) . So, the above definition is a generalisation of Definition 13. II 

Example 3.3 (Guarded composition of measurements) We consider two simplest measure- 
ments; that is, measurements on a qubit in the computational basis |0), |1) and in basis 
|±> = ^(|0>±|1»: 

M° = {M ° = |0)(0|,Af 1 ° = |l)(l|}, 
M l = {Ml = |+)(+|, M\ = |-}(-|}. 
Then their guarded composition along another qubit is measurement 

M = (\0)^M )D(\l)^M 1 ) 
= {M 00 ,M 01 ,M W ,M 11 } 

on two qubits, where 

M^(|^o)|0) + = i=(MP|^o)|0) + Mj\^)\l)) 

for any states \ipo), \ipi) of a qubit and i, j 6 {0, 1}. 

The following lemma shows that the guarded composition of operator-valued functions 
is well-defined. 

Lemma 3.2 The guarded composition F = D™ =1 \i) — > Fi is an operator-valued function 
in %(37-l s over ®" =1 Ej. In particular, if all Fj (1 < i < n) are full, then so is F. 
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3.4 Guarded composition of super-operators 

Guarded composition of a family of super-operators can be defined through the guarded 
compsition of the operator-valued functions generated from them. 

Definition 3.5 For each 1 < i < n, let £i be a super-operator in Hilbert space %. Let 7i s 
be a Hilbert space with {\i)} as an orthonormal basis. Then the guarded composition of £i 
(1 < i < n) is defined to be the family of super-operators: 

□? =1 \i) ->• Si = {£ (Df =1 \i) -> Fi) : 

Fi e ¥(£i) for every 1 < i < n}. 

It is easy to see that if n = 1 then the above guarded composition of super-operators 
consists of only S\. For n > 1, however, it is not a singleton, as shown by the following: 

Example 3.4 Let £q and £\ be the super-operators in Hilbert space % defined by unitary 
operators Uq, U\, respectively; that is, £; L = U\ o Uj (i = 0, 1). We set U to be the guarded 
composition ofUo and U±: U = (|0) — > L r o)D(|l) — > U\). Then the super-operator defined 
by U is £{U) G (|0) £oP(|l) -> Si)- Indeed, we have: 

(|0) -> £oM\l) -)• fx) = {£* = C/ e o Ul : < 9 < 2tt}, 

where Uq = (|0) — > f7o)D(|l) — > e %e U\). The non-uniqueness of the members of the above 
guarded composition is caused by the relative phase 9 between Uq and U\. 

4 Semantics of QGCL 

We first introduce several notations needed in this section. Let % and %' be two Hilbert 
spaces, and let E be an operateor in %. Then the cylindrical extension of E in % (g> %' is 
defined to be the operator E (8) I^i, where I^i is the identity operator in %' . For simplicity, 
we will write E for E ® whenever confusion does not happen. Let F be an operator- 
valued function in Ti over S. Then the cylindrical extension of F in % (g> "H' is the operator- 
valued function F' in % (8> over A defined by F'(5) = F(S) I-n> for every 5 6 A. For 
simplicity, we often write F for F' whenever confusion can be excluded from the context. 
Furthermore, let £ = J2i E% ° be a super-operator in H. Then the cylindrical extension 
of £ in H ® %' is defined to be the super-operator: £ = Ylii^i ® ^w) ° "X 1 ^H')- Here, 
for simplicity, the same symbol £ is used to denote the extension of £. In particular, if E 
is an operator in %, and p is a partial density operator in % <8> %' , then FpF^ should be 
understood as (F (g> i"^/)p(F^ (g> If £i and are two super-operators in a Hilbert 

space %, then their (sequential) composeition £\;£2 is the super-operator in H defined by 
(£i;<?2)(p) = £2{£\{p)) for any partial density operator p in %. 

4.1 Classical states 

We now define the states of classical variables in QGCL. 
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Definition 4.1 The (partial) classical states and their domains are inductively defined as 
follows: 



1. e is a classical state, called the empty state, and dom(e) = 0; 

2.1fx£ Var is a classical variable, and a G D x is an element of the domain of x, then 
[x <— a] is a classical state, and dom([x ^— a]) = {x}; 

3. If both 5\ and 82 are classical states, and dom(8\) n dom(82) = 0, then 8182 is a 
classical state, and dom(8i82) = dom(8\) U dom(82); 

4. If Si is a classical state for every 1 < i < n, then ©™ =1 <5j is a classical state, and 

n 

dom(®f =1 8i) = [J dom(8i). 
i=i 

Intuitively, a classical state 8 defined by clauses (1) to (3) in the above definition can 
be seen as a (partial) assignment to classical variables; more precisely, 8 is an element of 
6 6 Uxedom(S) Dx', that is, a choice function: 8 : V -> {J x&dom{S) D x such that 8{x) G D x 
for every x G dom(8). In particular, e is the empty function. Since Y\ xe a D x = {e}, e is 
the only possible state of with empty domain. The state [x a] assigns value a to variable 
x but the values of the other variables are undefined. The composed state 5\ 82 can be seen 
as the assignment to variables in dom{8{) U dom(82) given by 

{8 1 82)(x) = { h{x) if * ed0 "fi)> (14) 
I £2(2;) if x G dom(82)- 

Eq. (fl4l) is well-defined since it is required that dom{8\) n dom(82) = 0- In particular, 
e5 = (5e = 8 for any state 8, and if x ^ dom(8) then 5[x ^— a] is the assignment to variables 
in dom(8) U {x} given by 



<5[x <- a](y) 



8(y) if y G dom(8), 
a if y = x. 



The state ®" =1 5i defined by clause (4) in Definition 14. 11 can be thought of as a kind of 
superposition of 5{ (1 < i < n). 



4.2 Semi-classical denotational semantics 

For each QGCL program P, we write A(P) for the set of all possible states of its clas- 
sical variables. The semi-classical denotational semantics \P~\ of P will be defined as an 
operator- valued function in ~H qvar (p) over A(P), where T-L qva r(P) i s tri e type of quantum 
variables occurring in P. In particular, if qvar(P) = 0; for example P = abort or skip, 
then Ti qvar ^p) is a one-dimensional space and an operateor in H.^ can be identified 
with a complex number; for instance the zero operator is number and the identity operator 
is number 1. For any set V C qVar of quantum variables, we write Iy for the identity 
operator in Hilbert space Hy 
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Definition 4.2 The classical state A(P) and semi-classical semantic function [P] of a 
QGCL program P are inductively defined as follows: 

1. A(abort) = {e}, and [abort] (e) = 0; 

2. A(skip) = {e}, and [skip] (e) = 1; 

3. A(U[q]) = {e}, and \U\qJ] (e) = Uq, where Uq is the unitary operator U acting in 

Hqi 

4. IfP = M[x <- q] : {P m }, where M = {M m }, then 

A(P) = \J{S[x ^m]:Se D(P m )}, 

m 

[P] (6[X <- m\) = ([P m ] (5) ® I V \q Var{Pm) ) ■ (M m ® Iy\g) 

for every 5 G A(P m ) and for every m, where V = q U |J m qvar(P m ); 

5. IfP = D? =1 q,\i)-tP if then 

n 

A(P) = 0A(P), 

i=i 

\P]=DU |i)^[P]; 

6. 

A(P i; P 2 ) = A(Px);A(P 2 ) 

= {StS 2 : 5i G A(Px) and <5 2 G A(P 2 )}, 

[Pi! P 2 ] = ([P 2 ] (5 2 ) 8) /v\ gW r(P 2 )) 

• ({PlWl) ® Iv\qvar{P!)) 

where V = qvar{P\) U qvar(P2); 

Since it is required in Definition 12. 1 1 that far (Pi) n war(P 2 ) = in the sequential 
composition Pi; P 2 , we have dom{8\)C\dom(82) = for any 5\ G A(Pi) and 5 2 G A(P 2 ). 
Thus, Eq. (fT5T > is well-defined. Intuitively, if a quantum program P does not contain any 
guarded command, then its semantic structure can be seen as a tree with its nodes labelled 
by basic commands and its edges by linear operators. This tree grows up from the root 
in the following way: if the current node is labelled by a unitary transformation U, then a 
single edge stems from the node and it is labelled by U; and if the current node is labelled 
by a measurement M = {M m }, then for each possible outcome m, an edge stems from 
the node and it is labelled by the measurement operator M m . Thus, each classical state 
6 G A(P) is corresponding to a branch in the semantic tree of P, and the value of semantic 
function [P] in state 5 is the (sequential) composition of the operators labelling the edges 
of 5. The semantic structure of a quantum program P with guarded commands is much 
more complicated. We can imagine it as a tree with superpositions of nodes that generate 
superpositions of branches. The value of semantic function [P] in a superpositions of 
branches is then defined as the guarded composition of the values in these branches. 
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4.3 Purely quantum denotational semantics 

Now the purely quantum semantics of a quantum program can be naturally defined as the 
super-operator induced by its semi-classical semantic function. 



Definition 4.3 For each QGCL program P, its purely quantum denotational semantics is 
the super-operator [P] in T-L qvar (p) defined as follows: 

lPj=£(\P])= £ \P](6)o\P](6)l 
5eA(P) 

The following proposition presents a representation of the purely quantum semantics of 
a program in terms of its subprograms. 



Proposition 4.1 1. [abort] = 0; 

2. [skip] = 1; 

3. [P 1 ;P 2 ] = [P 1 ];[P 2 ]; 

4. {U[q]\ = U 1 oU 1 ; 

5. \M[x «- q] : {P m }] = £ m [(M m o Mm); [P m ]] . Here, [P m ] s/iokW fee seen as a 

cylindrical extension in T-Ly from T~Lq V ar(P m )> M m o M m as a cylindrical extension in 
%y from Hq, and V = q U \J m qvar(P m ); 



6. Qy K) — ^ -Pil £ 0™ =1 |i) — )■ [PjJ.^ere [Pj] should be understood as a cylindri- 

cal extension in T-Ly from ~rlq V ,ar{P l )f or every 1 < i < n, and V = (?U|J™ =1 qvar(Pi). 



The symbol "G" in clause 6) of the above proposition can be understood as a refinement 
relation. It is worth noting that in general "e" cannot be replaced by equality. This is exactly 
the reason that the purely quantum semantics of a program has to be derived through its 
semi-classical semantics and cannot be defined directly in a compositional way. 

Equivalence relation between quantum programs can be introduced based on their purely 
quantum semantics. 

Definition 4.4 Let P and Q be two QGCL programs. If qvar(P) = qvar(Q) and [P] = 
[QJ, then we say that P and Q are equivalent and write P = Q. 

4.4 Weakest Precondition Semantics 

The notion of quantum weakest precondition was proposed by D'Hondt and Panangaden |5). 



Definition 4.5 Let P be a program, and let M and N be positive (Hermitian) operators in 

Hqvar(P)- 
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1. Iftr(Mp) < ir(iV[P] (p))for all p G T^{T~L q var(P))' then M is called a precondition 
of N with respect to P. 

2. N is called the weakest precondition of M with respect to P, written N = wp.P.M 
if 

(a) N is a precondition of M with respect to P; and 

(b) N' C N whenever N' is a also precondition of M with respect to P. 

wp.P can be seen as the super-operator in T-L qva r(P) defined as follows: for any positive 
operator M, (wp.P)(M) = wp.P.M is given by clause 2) of the above definition, and 
wp.P can be extended to the whole space C{T-L qvar ^) by linearity. 

The weakest precondition semantics of QGCL programs are given in the next proposi- 
tion. 

Proposition 4.2 For any QGCL program P, and for any positive (Hermitian) operator M 
in l-L qV ar(P)> wp.P.M is given as follows 

1. wp. abort = 0; 

2. wp.sk.ip = 1; 

3. wp.(Pi;P2)=wp.P 2 ;wp.Pi; 

4. wp.U[q] = U\ o U q ; 

5. w P .(M[x<-q]:{P m }) = Y m 

6. wp.{0? =l q, \i) Pi) e Df =1 \i) -> wp.Pi. 

Some cylindrical extensions of super-operators are used but unspecified in the above 
proposition because they can be recognised from the context. Again, "G" in the above 
clause 6) cannot be replaced by equality. 

5 Quantum Choices: Superpositions of Programs 
5.1 Definition and Example 

As explained in Sec. [T] quantum choice may be defined based on quantum guarded com- 
mand. 

Definition 5.1 Let P and Pi be programs for all 1 < i < n such that q = qvar(P). If 
{\i}} is an orthonormal basis ofH q , and q n UiLi QVar(Pi) = 0, then the quantum choice 
of Pi, P n according to P along {\i)} is defined as 

n 

P,\i) ^ Pi = P;n? =1 q,\i) ^ Pi- 

i=l 



wp.P m ; (Mjn O Mm) 
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In particular, if n = 2, then the quantum choice will be abbreviated to P\ prgi © P 2 or 
Pip® P 2 . 



Example 5.1 Quantum walks have been extended to include multiple walkers and coins. 
These extended quantum walks can be conveniently written as QGCL programs with quan- 
tum choice. We consider two quantum walkers on a line sharing coins H16\l . The Hilbert 
space of a single walker on a line is H = Tip <8> %& where H p = span{\x) : x £ 
Z (integers)} is the position space and 7i c = span{L, R} is the coin space. Its step 
operator is W = (Tl <S> \L)(L\ + Tr (g> \R)(R\)(Iu (g) H), where 1% is the identity 
operator in Hp, 



is the 2x2 Hadamard matrix, and Tl , Tr are left- and right-translation, respectively; that 
is, Tl\x) = \x— 1), Tr\x) = \x+ 1) for every x £ Z. Then the Hilbert space of two walkers 
is H % 7~L, and if the two walkers are independent, then the step operator is W ® W. A 
two-qubit unitary operator U can be introduced to entangle the two coins and it can be 
thought as that the two walkers are sharing coins. A step of two walkers sharing coins can 
be written as a QGCL program as follows: 



where q± , q 2 are the position variables and c\ , c 2 the coin variables of the two walkers, 
respectively. 

5.2 Local Quantum Variables 

A quantum choice is defined as a "coin" program followed by a quantum guarded command. 
A natural question would be: is it possible to move the "coin" program to the end of a 
guarded command? To answer this question positively, we need to extend the syntax of 
QGCL by introducing block command with local quantum variables. 

Definition 5.2 Let P be a QGCL program, let q C qvar(P) be a sequence of quantum 
variables, and let p be a density operator in Hq. Then 

I. The block command defined by P restricted toq = p is: 



3. The purely quantum denotational semantics of the block command is give as follows: 




U[c 1 ,c 2 }; {T L [qi}H[ Cl ] © T R [qi\); (T L [q 2 ] H[c2] © T R [q 2 ]) 



begin local q := p;P end. 



2. The quantum variables of the block command are: 



qvar (begin local q := p;P end) = qvar(P) \ q. 



[begin local q := p; P end] (a) = tr« ? ([P](a © p)) 



for any density operator a in H, 



■qvar(P)\q- 
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The following theorem shows that the "coin" in a quantum choice can be move to the 
end of the guarded command if encapsulation in a block with local variables is allowed. 



Theorem 5.1 

n 

Qu[q], \i) ^Pi^ (qtif^K) -»• Pi)-,u[i]- (i6) 

8=1 

More generally, we have: 

n 

ffi-P, \i) ->■ Pi = begin local r := \ip ); 

)ti (17) 
□ijlV'ii) -> <9ii! end 

for some new quantum variables r, state \(fo) € orthononnal basis {\ipij)} ofHq®Hr, 
programs Qij, and unitary operator U in Tig <£> Hr, where q = qvar(P). 



5.3 Quantum implementation of probabilistic choices 



We now examine the relation between probabilistic choice and quantum choice. To this end, 
we first extend the syntax of QGCL by adding probabilistic choice. 

Definition 5.3 Let Pi be a QGCL program for each 1 < i < n, and let {pi}" =1 be a 
sub-probability distribution; that is, p$ > for each 1 < i < n and Y17=iPi — 1~ Then 

1. The probabilistic choice ofP\,...,P n according to {pi}" =1 is 

n 

^Pi@Pi- 
i=l 

2. The quantum variables of the choice are: 

/ n \ n 

qvar y~] Pj@Pi = qvar(Pi). 



i=i 



3. The purely quantum denotational semantics of the choice is: 



i=l 



'-Pi 



E 

i=i 



Pi 



Example 5.2 (Continuation o f Example \3.3\ Probabilistic mixture of measurements) It is 
often required in quantum cryptographic protocols like BB84 to randomly choose between 
the measurement M° on a qubit in the computational basis and the measurement M 1 in 
the basis |±). If we perform measurement M' 1 on qubit and discard the outcomes 
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of measurement, then we get pi = Mq\^){^\Mq + M{\ip)(ip\M{ for i = 0, 1. We now 
consider the unitary matrix 

on a qubit, where p, q > and p + q = 1. Let 

P = begin local q : = |0); g := 

□i=o,i9, K) -> end 

where q, q\ are qubit variables. Then for any \ip) 6 a«J i, j G {0, 1}, we Ziave: 

|Y^) £ Af y (|^>tr|o» = yf m<V)|o> + y|M/|^)|i), 
[P](|V>M) = tr*, I ^ 

\i,3=0,l J 

«J=0,1 
= PPo + 9Pi- 

5(3, program P can be seen as a probabilistic mixture of measurements M° and M 1 . 

As shown by the following theorem, if the "coin" variables are treated as local variables, 
then a quantum choice degenerates to a probabilistic choice. 

Theorem 5.2 Let qvar(P) = q. Then we have: 

n n 

begin local q := p; P, \i) — > Pi end = ^ Pi@p% 
i=i i=i 

where pi = (i\\P\(p)\i) for every 1 < i < n. 

Conversely, for any probability distribution {pi}f =1 , we can find annxn unitary op- 
erator U such that pi = \Uio\ 2 (1 < i < n). So, it follows immediately from the above 
theorem that a probabilistic choice J27=i Pi®Pi can always be implemented by a quantum 
choice: 

n 

begin local q := |0); ^ U[q], \i) ->■ Pi end 

i=i 

where q is a family of new quantum variables with an n— dimensional state space. 
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6 Quantum Recursion 



Now we need to further extend the syntax of QGCL. We first add a countable set of program 
names, ranged over by X, Y, to the alphabet of QGCL, and then introduce the following: 

Definition 6.1 QGCL programs are defined by combining Definitions \2.1\ \5?2\ \5.3\ and the 

following two clauses: 

1. Every program name X is a program, and both var(X) and qvar(X) are given a 
priori. 

2. IfP is a program and X a program name such that var{P) C var{X) and qvar(P) C 
qvar(X), then pX.P is a program, and var(fiX.P) = var(X), qvar(fiX.P) = 
qvar(X). 

We consider a special case of quantum recursion, namely quantum loop, and show an 
interesting difference between quantum loops with classical control flows defined in lfl8l 
and quantum loops with quantum control flows. The quantum loops considered in lfl"8l can 
be written as QGCL programs of the form: 

Loop = while M[q] = 1 do q := Uq 

= pX.M[x <- q] : {P = skip, P x = q := Uq; X} 

where q is a sequence of quantum variables, M = {Mo, Mi} a binary ("yes-no") measure- 
ment in T-lq and U a unitary operator in 7iq. The control flow of Loop is determined by 
measurement M in the loop guard: if the outcome of measurement is then Pq is executed 
- the loop terminates; if the outcome of measurement is 1 then Pi is executed - the program 
executes the loop body q := Uq and then runs into the loop again. Program Loop can be 
approximated by a series of iterations {Q n }^ =0 defined as follows: 

{Qo = abort, 
Q n+1 = M[x <- q] : {P = skip, (18) 
= q : = Uq; Q n } (n > 0). 

If the classical control flows of Q n (n > 0) determined by the outcomes of measurement 
M are replaced by quantum control flows defineed by quantum choices, then we obtain the 
following quantum iterations: 

{Q' Q = abort, 
Q'n+l = ski Pc[ Qn+1 ] © (V ■= Uq; Q' n ) (n > 0) 

where C is a "coin" 2x2 unitary matrix. It is worth noting that we have to introduce a 
sequence q\ , q2 , . . . of new qubit variables in order to well-define the quantum choices used 
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in Q' n (n > 1). For each n > 0, since var{Q' n ) = and qvar(Q' n ) = q U {q%, q n }, 
the semi-classical semantics \Q' n ~\ of Q' n is an operator- valued function in T~Lq Cg) ®" =1 
over {e}. Suppose that the input is a state \ip) in and all the auxiliary qubit variables 
gi , . . . , q n are initialised in state 1 0) . For simplicity of calculation, we take C = H (the 2x2 
Hadamard matrix; see Example 15. lb . Then 

n— 1 1 

rQ' n i(6)^)ior = ^-==c/^)|o)™-|ir. 

It is clear that we cannot directly define the semantics of a quantum loop as the limit of 
{Q' n }^=o because the state spaces of Q' n are different for different n. To overcome this 
difficulty, a natural idea is to localise qubit variables q\ , . . . , q n : 

Qn = begin local q u q n := \0) n ; Q' n end. 

But such a localisation makes the quantum iterations degenerate to probabilistic iterations: 

n— 1 -. 

i=0 

This gives an example showing that quantum loops, or more generally quantum recursions, 
with quantum control flows are much harder to deal with than those with classical control 
flows. Due to the limited space, a more detailed treatment of quantum recursion is post- 
poned to another paper. 

7 Conclusions 

Three new quantum program constructs - quantum guarded command, quantum choice and 
quantum recursion - are defined in this paper. We believe that introducing these constructs 
is a significant step toward the full realisation of "quantum control" in quantum program- 
ming. In the further studies, we will consider quantum recursions with quantum controls in 
detail, and we will establish various algebraic laws for QGCL programs that can be used in 
program transformations and compilation. A quantum Floyd-Hoare logic was built in ifTTl 
for quantum programs with only classical control flows. So, another interesting topic for 
further studies would be to extend this logic so that it can also be applied to programs with 
quantum control flows. 
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Appendix: Proofs 

7.1 Proof of Lemma 03| 

We start with an auxiliary equality. Put: 

(5ieAi,...,5„eA n 
For any \<p), G % <g> H s , we can write: 

n 

i=l 
n 

w = Z>>i<> 

i=l 

where |^), |^) G % for each 1 < i < n. Then 

Si,...,S n 

= e i: (n*.) (n*»« 

<5i,...,<5„ i,i'=l J \k^i' 



{iPiWiiStfFviSi,)^)^') 




=t e 

i=l <5i,...,5i_i,5 i+ i,...,(5 n \kj^i 

k 

n 

= EE^i F ^) tF ^)i^> 

i=l Si 
n 

= E^iE i? ^) ti? ^)i^> 

i=l <5i 

because for each A;, we have: 

E i A **J 2 = lj 

and thus 
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(1) We now prove that F is a semi-classical semantic function in 'H®'H S over ©™ =1 A n . 
It suffices to show that F C In®~H 3 '' that is, (^|F|</?) < for each \tp) £ 7~L (3 7i s . In 
fact, for each 1 < i < n, since F{ is a semi-classical semantic function, we have: 

Si 

{<Pi\^2 Fiidrf Fi(5i)\ipi) < (tpi\(pi). 
Then it follows from Eq. (fT9l ) that 

n 

{<P\F\<P) < J^^ilVi) = (vlv)- 

i=l 

(2) For the case where all Fi (1 < i < n) are full, we have: 

Si 

for all 1 < i < n, and it follows from Eq. ( fl9l ) that 



i=i 



So, it holds that F = In®n s by arbitrariness of |</?) and -0), and F is full. 

7.2 Proof of Proposition O 

Clauses 1) - 4) are obvious. 

5) By definition, for any partial density operator p in V, qvar {P)^ we have: 

{M[xq]:{P m }\{p) 
m <5eA(P m ) 

= ^ E (\Pm~\($) ® Iqvar(P)\qvar(P m )) 
m <5eA(P m ) 

(M m <g> I gV ar{P)\q)p( M L ® Iqvar(P)\q) 
(rP m l(<5) f (8) Iy\g mr (p m )) 

=E E (r^miw® I qvar(P)\qvar(P m ) ) 
m <5eA(P m ) 

(M m/ oA4)(rP m l(*)t ® I n , mr (P m) ) 

= ^[P m I(M m pA4) 
m 

= ^oMl);^] (p). 
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6) For simplicity of the presentation, we write: 



A 



By definition, we obtain: 



Since 



P = D? =1 q, \i) -)• Pi. 



lPj=S(\0^ =1 q,\i)^P l ]). 



and \Pi] £ F([Pi]) for every 1 < i < n,it holds that 

[Pfle |i}^P):PeF([Pfl) 

for every 1 < i < n} = E\™ =1 \i) - 



7.3 Proof of Proposition g3] 

The proof is based on the following key lemma by D'Hondt and Panangaden [5 ]. 

Lemma 7.1 If the semantic function [P] of program P has the Kraus operator-sum repre- 
sentation: [P] = ^ ■ Pj o Ej, then we have: wp.P = Y^j Ej ^j- 

Now we start to prove Proposition 14.21 Clauses 1) - 4) are immediate from Proposi- 
tion [4j] and Lemma 17711 

5) Suppose that for every m, 

m 

Then by Proposition 14. 11 5) we have: 

[M[x <- 5] : {P m }fl = £ [i M m o Ml); [P m j 



£ (M m oMt l );^(p mim oP ? t 
mi. i m 

m i m 
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Using Lemma 17711 we obtain: 

wp.(M[x <- q]; {P m }) = ^ E i E mi m M m )^ o (E mim M m ) 

m i m 



E 



[w P .P m ;(MloM n 



6) For each 1 < i < n, assume that the semi-classical semantics of Pi is the function 
\Pi] over A = {ji} such that 

\P i -\{j i ) = E ik 
for every ji. Then by Definition [43] we obtain: 

lP i l=Y,E iji oEl ji , 

ji 

and it follows from Lemma 177X1 that 

For any \cp) = Yh=x l<Pi)N)> where \ G ^-gvar[Pi) 0- < i < n )> we define: 



\ T,i k ( E li k ) iE li k \Ei k E l k E ki k 



trEl jk E kjk 



and Xkjf.'s are defined by Eq. (fT3T ). By Definitions I3.4l and l3.5l we have: 



On the other hand, by Definitions 14.21 (5) and !4.3l we have: 

[□2=1 9, \i) -> fil = E ^-in^l-i, 



(21) 
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where i^v.^'s are defined by Eq. (fl2"l ). Applying Lemma ITT1 once again, we obtain: 



■3l—3n- 

3l,—,3n 



So, we now only need to prove that 

for all jx, ...,j n . In fact, for any \tp) = Ya=i IVOK) with |^) G U qv ar(Pi) (1 < * < ")> it 
holds that 



i=l i=l 

= X)Ci(ivi>»% i i^» 

i 



n 



i=l i=l 

(\<P),Fh-jM) 



because £j's are real numbers, and it follows from Eq. (1211 ) that 

Ci = n Afcj fc 

for each 1 < i < n. Thus, we complete the proof. 
7.4 Proof of Theorem EE 

We first prove Eq. (fT6l) . Assume that [\Pj] is the operator-valued function over A$ such that 
\P i ] (<^) = F iSi for each <5j G Aj (1 < z < n). We write: 

P = D2 =1 c4|i)->P i . 
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Then for any |</>) = £™ =1 |^)K>, where \^) e H v (1 < i < n), and V = [f i=l qvar(Pi), 
we have: 



= \P](®« =1 6 i ) 
= \p-\(®ti5i) 



8=1 



rt / n 



3=1 Wi 



si=l 



Let Lif S and RHS stand for the left and right hand side of Eq. ( fT6l ), respectively. Then it 
holds that 

= E(n A ^J F ^(E^i^))^ 

j= i \k^j j \i=i ) 



n I n 



E Ew w 

j=i \i=i ) 

n n 

El^>E(c^li» 

i=l j=l 



rntjI^-^fiKe^OlEl^W)) 



K« = l 



Consequently, it follows that = and we complete the proof of Eq. (TT6l ). 

Now we are ready to prove Eq. (fTTT ). Since [P] is a super-operator in "H^-, there must be 
a family of quantum variables r, a pure state \tpo) G a unitary operator U in 7^ ® 
and a projection operator K onto some closed subspace /C of %r such that 

\P\{p) = tr nT (KU(p \vo)(vo\)U j K) 

for all density operators p in %g (see the system-environment model of super-operators, 
Eq. (8.38) in iTTTlD . We choose an orthonormal basis of /C and then extend it to an orthonor- 
mal basis {\j)} of H^- Define pure states = U^\ij) for all i,j and programs 



Qi 



[Pi if I J) G K, 

[abort if \j) <£ K. 
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Then by a routine calculation we have: 

|y) -> Q -](<r) = [□<!») -> PiK^aif) (22) 

for any <r G Hquruv> where 1/ = UILi Qvar(Pi). We now write i?i/S for the right hand 
side of Eq. (fTTT ). Combining Eqs (fT6l ) and (l22l . we obtain: 

[p#s](p) 

= tr %Y ({UijU^ij) -»■ Qij-;?7[g,r]](p® l^oXvo 



= ^ [©^M.ltf) -> Qyl(p® |<Po)<<Po|) 

= trur (Pylu) ->■ l^oXvol 
= tr^pili) -> P](itt/(p® |<p )M)^) 
= [Dili) -> P^r^ifC/^® l^oXvol)^^)) 

= [□<!») -+Pi]([P](p)) 
= [0p,K>^^1(p) 

i 

for all density operators p in 7{^. Therefore, Eq. (fTTT ) is proved. 

7.5 Proof of Theorem [53] 

To simplify the presentation, we write: 

and assume that [~Pj] is the operator- valued function over Aj such that \P{] (Si) = E i s i for 
each <5j G Aj. Let G gi ;ar(Pi) an d |<p) G 'Hg. We can write: 



n 



1=1 



for some complex numbers a% (1 < i < n). Then for any <5j G Aj (1 < z < n), we have: 

i^,.. a „) = rpi(®r=i5i)(i^)) 
= r«i(er=i^) f^oii^A 

= ^ ai n Afc ^ ^ 
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where Aj^ 's are defined as in Eq. ( fT3l >, 



\^ Sl ...s n )(^6 1 ...sJ = ^^a* ]Jhs k Hhs k 

i,j=l \k^i J \k^j 

E i5i \^){^\E] s .®\i){jl 



and it follows that 



i=l \k^i 



Using Eq. (I20I ). we obtain: 

tr n JR\(\^){ipM) 



\5i,...,8 n 

tr H -\^Su..s n )(^s l ...s n \ 

8i,...,8 n 



E 



a; 



E in** 



<5; 



8=1 



Now we do spectral decomposition for [P] (p) and assume that 

[pj(p) = 5>h><^i- 



We further write: 



m) 



(23) 



i 

for every Z. For any density operator a in 'H(j n _ 1 g t ,ar(P i )» we can write a in the form of 

<y = y~y m |^ m )(^ m |. 
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Then using Eq. (1231 . we get: 



where 



[begin local q := p; U^ =l P, \i) -)■ Pi end]](<7) 

= tr H J{P;R\(a®p) 
= tr H JRj(a®lP}(p)) 



tr-HqlRj ^r m si\ip m ipi){ipiip m \ 

\ m,l 

E r mSitr ni {Rj (\ip m ipi) (<piip m \) 

rnj 

n 

^r mS ^|a K | 2 [Pi](|^ m )(^ m |) 



J i=l 



l i=l \ m 

n 

I i=l 



i=l \ I 

n 



i=l 



Pi = ^silaul 2 = y^si(i\(pi){ipi\i) 
i i 

= (i\ fe-ilwKwl) \ i ) = (*\l p l(p)\ i )- 
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